Energy size effects of two-dimensional Ising spin glasses 



I. A. Campbell, "'^ Alexander K. Hartmann,^ and Helmut G. Katzgraber'^ 

' Laboratoire des Verres, Universite Montpellier II, 34095 Montpelher, France 
^Institut fiir Theoretische Physik, Universitat Gottingen, Tammannstrafie 1, 37077 Gottingen, Germany 
Theoretische Physik, ETH Honggerberg, CH-8093 Zurich, Switzerland 
(Dated: February 6, 2008) 

We analyze exact ground-state energies of two-dimensional Ising spin glasses with either Gaussian 
or bimodal nearest-neighbor interactions for large system sizes and for three types of boundary 
conditions: free on both axes, periodic on both axes, and free on one axis and periodic on the 
other. We find accurate values for bulk-, edge-, and corner-site energies. Fits for the system with 
Gaussian bonds are excellent for all types of boundary conditions over the whole range of system 
sizes. In particular, the leading behavior for nonfree boundary conditions is governed by the stiffness 
exponent 9 « —0.282 describing the scaling of domain-wall and droplet excitations. For the system 
with a bimodal distribution of bonds the fit is good for free boundary conditions but worse for other 
geometries, particularly for periodic-free boundary conditions where there appear to be unorthodox 
corrections to scaling up to large sizes. Finally, by introducing hard bonds we test explicitly for the 
Gaussian case the relationship between domain walls and the standard scaling behavior. 

PACS numbers: 75.50.Lk, 75.40.Mg, 05.50.-fq 



I. INTRODUCTION 

The properties of the two-dimensional Edwards- 
Anderson Ising spin glass (ISG)i^iMi^ with either Gaus- 
sian or bimodal (± J) nearest-neighbor interactions have 
been intensively studied. Efficient algorithmsSiS have 
been developed which provide exact ground-state ener- 
gies and configurations for these systems up to large sys- 
tem sizes with N — LxL spins up tc^ L — 1800, allowing 
for accurate measurements of fundamental parameters 
and stringent tests of theoretical predictions including 
precise evaluations of size-effect corrections to scaling. 

In both Gaussian and bimodal disorder distributions, 
freezing occurs only at zero temperaturefiSi i.e. Tc — 0. 
For the Gaussian ISG as for any ISG system with a 
continuous distribution of interactions there is a unique 
ground state (together with its reversed image). At 
T = the thermodynamic limit correlation function G{r) 
is constant, G(r) = 1 for all distances r. As the critical 
exponent 7] is defined through G(r) ^ 7.-(<i-2+»7) at the 
ordering temperature, ry = 0. For the ± J ISG, the ground 
state is highly degenerate, with a residual entropy at zero 
temperatureii S{T = 0) = 0.078(5)fcB. In consequence 
the time averaged correlation function G(r) drops with 
increasing r, meaning rj must be positive^ The domain- 
wall stiffness has been measured accurately through com- 
parisons of the energies of samples with periodic and 
antiperiodic boundary conditions in one direction. For 
the Gaussian ISG, 6 = -0.282(2) From the scal- 
ing rules applied at a zero temperature transition^^ii^ 
we obtain for the critical exponent ly = —I/O, hence v = 
3.55(3). 1^-^^ For the ± J 2d ISG, large-L simulations^o on 
samples with periodic/antiperiodic boundary conditions 
along one axis and free boundary conditions on the other 
axis showed 9 — 0, with significant corrections to scaling 
up to L ~ 100. Finite-size scaling and Migdal-Kadanoff 
simulations on a variety of 2d ISG samples with different 



distributions of interactions arc consistent with a unique 
lower critical dimension di ^ 2.5; for all the samples 
with discrete distributions of interactions 6* = at all di- 
mensions below the lower critical dimension^ For these 
systems the effective exponent v is infinite. The cor- 
relation length diverges exponentially^^ and the specific 
heat drops to zero as exp(— n/3); data were analyzed us- 
ing n equal tc^-^J^ 4 but recent work^" is consistent with 
n — 2. Although we will be concerned here only with 
zero-temperature behavior, we can note that there have 
been conflicting estimates of the "droplet" excitation ex- 
ponent for the Gaussian ISG, which standard arguments 
indicate should be identical to 9. It has been suggested 
that apparent disagreements are artifacts due to correc- 
tions to scalingSi but Berthier and Young foundSS that 
whether such corrections are visible seems to depend on 
the details of the measurements. 

We have carried out measurements of the ground-state 
energy per spin on square samples as a function of sys- 
tem size L for three different types of boundary condi- 
tion: free along both x and y pairs of boundaries (to be 
referred to as ffbc) , free along one pair of boundaries and 
periodic along the other pair to give a cylindrical geome- 
try (pfbc), and periodic along both pairs of boundaries to 
give a toroidal geometry (ppbc). These three geometries 
are physically distinct. In the ffbc case, there are bulk 
sites, edge sites, and corner sites. However the system 
has no boundary condition constraints, or in other words 
it can be taken to contain no "domain walls," at least 
for a sample with a unique ground state. In the pfbc 
scenario there are edge sites but no corner sites; on the 
other hand the sample is constrained by the boundary 
condition along one direction. Finally for the familiar 
ppbc case there are no edges or corners but there are 
constraints along both directions. It is for this geometry 
with no boundaries that the standard scaling approach 
should apply. 



2 



The behavior of the fFbc systems can be explained in 
simple geometrical terms. The behavior of the Gaussian 
systems with pfbc and ppbc is governed by an algebraic 
term with the exponent 9 « —0.282, the same value as 
the exponent describing the behavior of domain- wall and 
droplet excitationsi2iS4 For this reason, we also com- 
pare sample by sample ffbc and pfbc realizations with 
the same interaction set and look for the appearance 
of system-spanning domain walls. Finally, we introduce 
hard bonds to force system-spanning domain walls and 
study how the appearance of these depends on the frac- 
tion of hard bonds. 

The paper is organized as follows: In Sec. ^ we in- 
troduce the model, observables, and details of the al- 
gorithms used. Scaling arguments are summarized in 
Sec. mil Results for the ffbc, ppbc, and pfbc scenar- 
ios are presented in Sees. IIVI Ivl and IVII respectively. 
Domain-wall calculations are discussed in Sec. IVIII and 
concluding remarks are contained in Sec. IVIIII 



II. MODEL AND ALGORITHMS 

The Hamiltonian of the two-dimensional Ising spin 
glass is given by 

H. = —"^^ JijSiSj, (1) 

where the sites i lie on a square lattice in two dimen- 
sions and the Jij are nearest-neighbor interactions. In 
the Gaussian case the couplings are chosen according to 
a Gaussian distribution with zero mean and standard de- 
viation unity, whereas in the bimodal case the couplings 
Jij can take values of {±1} with equal probability. 

In greater than two dimensions, or in the presence of a 
magnetic field, the exact calculation of spin-glass ground 
states belongs to the class of A^P-hard problems 
This means that only algorithms with exponentially in- 
creasing running times are known. Here, we have studied 
mainly square lattices with periodic boundary conditions 
in at most one direction without external magnetic fields. 
For this special case of a planar system there are effi- 
cient polynomial-time "matching" algorithmsi^ The ba- 
sic idea is to represent each realization of the disorder by 
its frustrated plaquettes,-'' Pairs of frustrated plaquettes 
are connected by paths in the lattice and the weight of 
a path is defined by the sum of the absolute values of 
the coupling constants which are crossed by the path. A 
ground state corresponds to the set of paths with mini- 
mum total weight, such that each frustrated plaquette is 
connected to exactly one other frustrated plaquette. This 
is called a minimum-weight perfect matching. The bonds 
which are crossed by paths connecting the frustrated pla- 
quettes are unsatisfied in the ground state, and all other 
bonds are satisfied. 

For the calculation of the minimum-weight per- 
fect matching, efficient polynomial-time algorithms are 



TABLE I: The number of samples for the ground-state cal- 
culations for different system sizes, distributions and free-free 
(ffbc), resp., free-periodic (pfbc) boundary conditions. 

L Gaussian ffbc Gaussian pfbc ±J ffbc ±J pfbc 

2 1 X lO'' 4 X 10^ 

3 1 X 10** 1 X 10'^ 1 X 10^ 1 X 10*^ 

4 9.1 X 10^ 1 X 10^ 1 X 10** 9.4 x lO'^ 
6 5.1 X lO'^ 9.9 X 10** 5 X 10^ 5.4 x 10^ 
8 2.1 X 10^ 9.9 X 10** 2 X 10^ 3.4 x 10^ 
10 1.1 X 10^ 9.9 X 10** 1 X 10^ 1.4 X lO'^ 
16 4 X 10'* 9.9 X 10* 1 X 10* 4 x 10* 
20 1 X 10** 1 X 10** 1 X 10* 4 X 10** 
30 1 X 10** 1 X 10** 1 X 10* 4 X lO'* 
40 1 X 10** 1 X 10** 1 X 10* 3.9 X 10** 
60 1 X 10** 1 X 10** 1 X 10* 4 X lO'* 
80 1 X 10** 1 X 10" 1 X 10* 4 X lO'* 
120 1 X 10** 1 X 10** 1 X 10* 4 X 10** 
160 1 X 10** 1.1 X 10* 1 X 10* 4 X lO'* 
240 724 2800 5 x 10^^ 4.6 x 10* 
320 3710 5 X 10^ 2.4 x 10** 
480 2 X lO'* 



availableii^ Recently, an implementation has been 
presented)^ where ground-state energies of large systems 
of size N < 1800^ were calculated. Here, an algorithm 
from the LEDA library^* has been applied, which limits 
the system sizes due to the restricted size of the main 
memory of the computers which we used. 

Furthermore, we have studied two-dimensional sys- 
tems with fully periodic boundary conditions. Here, we 
use the spin-glass ground-state server at the University 
of Cologne^S where a "branch-and-cut" method22i2i (see 
Ref. for a tutorial on optimization problems and tech- 
niques, including branch-and-bound and branch-and-cut) 
is used, which is currently the fastest exact algorithm 
for computing spin-glass ground states;^ with the ex- 
ception of the polynomial-time special cases mentioned 
above. Nevertheless, the CPU time increases faster than 
a power of the system size. Therefore we cannot study 
as large systems as are possible with the matching algo- 
rithm. Nonetheless, the implementation of the branch- 
and-cut algorithm on the Cologne server is very efficient, 
hence we can still study quite a large range of sizes, in 
practice L < 64. 

The numbers of samples over which averages are taken 
are shown in Tables U and ^ 



III. SCALING 

We want to find the size-dependence of the energy from 
simple scaling arguments for the generic ppbc case and 
follow closely along the lines of the scaling discussion pre- 
sented in Ref. The internal energy per spin e of a 



1 X lO'' 




4 


X 


10® 




1 X lO" 


1 X 10 


1 


X 


lO'' 


1 X 10 


9.1 X 10^ 


1 X 10^ 


1 


X 


10® 


9.4 X 10^ 


5.1 X 10^ 


9.9 X 10* 


5 


X 


10^ 


5.4 X 10^ 


2.1 X 10 


9.9 X 10* 


2 


X 


10^ 


3.4 X 10 


1.1 X 10^ 


9.9 X 10* 


1 


X 


10^ 


1.4 X 10^ 


4 X 10 


9.9 X 10* 


1 


X 


10 


4 X 10* 


1 X 10** 


1 X 10* 


1 


X 


10* 


4 X 10* 


1 X 10** 


1 X 10* 


1 


X 


10* 


4 X 10* 


1 X 10** 


1 X 10* 


1 


X 


10* 


3.9 X 10* 


1 X 10** 


1 X 10* 


1 


X 


10* 


4 X 10* 


1 X 10** 


1 X 10* 


1 


X 


10* 


4 X 10* 


1 X 10** 


1 X 10* 


1 


X 


10* 


4 X 10* 


1 X 10** 


1.1 X 10* 


1 


X 


10* 


4 X 10* 


724 


2800 


5 


X 


10^5 


4.6 X 10* 




3710 


5 


X 


10^ 


2.4 X 10* 
2 X 10* 



3 



TABLE II: The number of samples for the ground-state cal- 
culations for full periodic (ppbc) boundary conditions, for dif- 
ferent system sizes and disorder distributions. 
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thermodynamic system can be expressed as a derivative 
of the free energy per spin / with respect to the temper- 
ature: 



1 d{Pf) 
'(3^ dT ' 



(2) 



where (3 — 1/T, T the temperature. 

We assume that the scaling of the temperature de- 
pendence of the correlation length is ^ ~ (T — Tc)'" 
or (T - Tc) - ^-^Z". Therefore 



-(1+1/^) 



Hence, we obtain 
e(r): 



1 djpf) 
/32 d^ dT 



B{T) 



(3) 



(4) 



B{T) is the smooth noncritical background. 

Using basic finite-size scaling arguments?^ we know 
that for the singular part of the free energy (3f ^ 

First we discuss the case Tc > 0. In this case the factor 
l//3^ is not critical. Therefore we obtain 



Pc 



(5) 



where /3c — l/Tc. At Tc, — CiL, Ci a constant, so the 
critical behavior of cl is 



Next, we discuss the case Tc 
P ~ ^1/"^, hence 



(Tc > 0) . (6) 
0. Now (3 is critical. 
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For a finite system of size L, 



B 



(Tc = 0) 



(7) 



(8) 



For Tc = one also has the scaling relation'^^ between 
the domain-wall stiffness exponent 6 and the exponent v, 
6 = — therefore 



-id-B) 



(Tc = 0) 



(9) 



For the ±J case, Tc = and the correlation length 
diverges exponentially. Following through the algebra 
witbiSi ^ = exp(2//3) leads to gl ~ Gcxj ^ L , which is 
just the same as Eq. © with 9 = 0, corresponding to 

= oo. 

Equation ijHJl is a scaling identity for any system with 
ppbc and Tc = 0, and there is no need to appeal to 
specific physical arguments to justify it. Any deviations 
from this identity must be due to corrections to scaling. 
(We can note that Bouchaud et alw^ refer to this scaling 
term as the "correction to scaling.") 

The leading correction to scaling should be either a 
renormalization group theory (RGT) irrelevant operator 
correction giving the form 



^(d-e) 



(1 +/CiL- 



(10) 



where uj is the RGT leading correction to scaling expo- 
nent and /Ci is a constant, or an analytic correction which 
plausibly introduces a term 1C2L~^ m^i^ This "analytic" 
term is unrelated to the RGT irrelevant operators. The 
arguments leading to a prediction for the energy cor- 
rection at this level can be very involved, even for the 
canonical 2c? Ising ferromagnet where there is no irrele- 
vant operator term but there are analytic terms in odd 
powers of 1/Li'^^i^^i'^° For ISGs the leading term in the 
RGT e-expansion for uj{d) as a function of dimension d 
isi^^: Lo{d) = (6 — c?) as we can safely assume that by 
d = 2, + w3>2 hence the dominating correction term 
is an analytic term. 

Thus finally, for ppbc from standard scaling arguments 
we expect 



(11) 



where the Ci term is the scaling term and the C2 term 
is the analytic correction term. In practice, if 9 is small 
compared to d, it will be hard to distinguish between this 
sum of two terms and an "effective" scaling. 



(12) 



with a single effective exponent 0off whose value is a func- 
tion of C2/ Ci and which can be larger or smaller than 6 
depending on the sign of this ratio. 

Although the remark is irrelevant for the two- 
dimensional case, we can note that if Tc > this dis- 
cussion is valid mutatis mutandi at Tc; however if Tc > 
scaling rules do not apply below Tc and in particular at 
T = 0. Alternative physical arguments must be used for 
discussing the ground-state behavior. 
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IV. FREE BOUNDARY CONDITIONS (FFBC) 

If boundary conditions are free along all four edges we 
can consider that there are no external constraints on the 
system. There are however edge and corner sites, edge 
bonds and corner bonds. The total number of bonds is 
equal to 2{L^ — L) while the number of spins is L^, hence 
we present a discussion in terms of the energy per spin 
CL or alternatively twice the energy per bond 2e^°"'* (the 
values will become identical in the infinite L limit). 

A heuristic geometrical scaling consists in writing twice 
the energy per bond as 

2ef'"^ = A*+B*/L + C*/L\ (13) 

where A*, B*, and C* are constants. Physically the first 
term in Eq. (|13l) represents the energy per bulk bond, the 
second term is a correction for the energy difference be- 
tween an edge and a bulk bond, and the third term a fur- 
ther correction for a corner bond. Even if the localization 
of the energy differences onto the edge or corner bonds 
is only approximate so the identification between terms 
and bonds is not rigorous, on geometrical grounds the 
ffbc energy size effects can be expected to have strictly 
this form with no further terms, at least down to small 
L values where neighboring corner effects begin to inter- 
fere. Translating Eq. H13|) into terms of energy per spin, 
i.e. using El = {L^ — L)e\°'^^ = L^e^, gives an energy 
per spin, 

BL^A*- {A* - B*)/L - {B* - C*)/L^ - C*/L^ . (14) 

For the isolated spin {L = 1) the energy per bond is not 
defined, but the energy per spin is identically zero. By 
inspection Eq. (|14|) happens to have a form such that cl 
must be exactly equal to zero for L = 1. We can also 
write 

eL = A + B/L + C/L^ +V/L\ (15) 

with A^ A*, B = -{A* - B*), C = -{B* - C*), and 
T> = -C*. For the Gaussian ffbc data the fit to Eq. ifT^ 
is shown in Fig.^ It can be seen that for the entire range 
of sizes from L = 1 to L = 240 the three-parameter fit 
is excellent with = 0.67. The same data with the 
same fit parameters A, B, C, and V are presented as a 
difference plot, i.e. as the difference between the fit and 
the actual data, in the inset of Fig. Remarkably, this 
naive equation with three free parameters is in agreement 
with the measurements to within one or two parts in 10^ 
for the whole range of L from L = 240 down to and 
including L = 1. The fit parameters are given in Table 

cni 

If we repeat the same analysis for ±J bonds, using 
data from L = 4 to L = 320, the agreement is almost 
equally good, cf. Fig. [3 The fit parameters (x^ = 1.29) 
are given in the Table UTTI (We have not included in the 
fit the point for L = 2 which lies almost 3% above the 
trend. By this size each system has only four bonds and 



TABLE III: Fit parameters according to Eq. lITll for ff 
boundary conditions, as well as for Gaussian and bimodal 
distributed disorder. 



parameter 


Gaussian 


±J 


A" 


-1.31479(2) 


-1.40197(2) 


B* 


-0.3205(9) 


-0.5492(20) 


C* 


0.042(3) 


0.506(18) 




1 IQi 102 

L 

FIG. 1: Energy per spin for the ffbc Gaussian case. The 
solid line represents a three-parameter fit to Eq. 11411 with 
= 0.67. The inset shows the difference between the fit and 
the actual data, Ae, represented as an energy per spin. Fit 
parameters are summarized in Table Hm 

there are only five possible energies for any ±J system. 
As one might foresee, quantization effects seem to break 
down the precise scaling by this size.) 

The fit extrapolates to a bulk (or infinite lattice size) 
energy per spin of 

e^^r - -1.31479(2) , (16) 

for the Gaussian case. The most recent published value 
is e^^^i'i^"'' = -1.317(1)30. For the ±J case the present 
data give a bulk energy per spin of 

Sbulk - -1-40197(2) . (17) 

This estimate is compatible with the very precise pub- 
lished value e^^^^ = -1.401938(2).^ Notice that the bulk 
energy of the bimodal system is significantly more nega- 
tive that the bulk energy of the Gaussian system. 

In the ffbc lattice with spins there are {L — 2)^ 
bulk sites, 4L — 8 edge sites, and 4 corner sites. We 
associate energies Cbuik, Cedgo, and ecornor with sites at 
each of these positions, respectively (equal to half the 
sum of the adjacent bond energies). Using El = L/^gl = 
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L 

FIG. 2: Energy per bond for the ffbc ±J case. The sohd 
hne corresponds to a three-parameter fit to Eq. 1141 with 
X'^ ~ 1.29. The inset shows the difference between the fit and 
the actual data. Fit parameters are summarized in Table ITTTl 



(i-2)^ebuik + (4L~8)eodgo+4ecornor and comparing with 
Eq. (|15|l . we can make the identifications for the energies 
per site, 



Cbulk ; 
eedgo = A+B/A, 
c-coincr 

A + B/2 + C/i. 



(18) 



Because each bulk site has four bonds, each edge site 
three bonds, and each corner site two bonds, and because 
each bond is shared with a neighbor site, to obtain the 
average bond energies for the bonds Unking to these three 
types of sites we divide Cbuik, Ccdgc, ecornor by 2, 1.5, and 
1, respectively. 

The bulk, edge, and corner energies per bond evalu- 
ated according to Eq. H18|) when inserting the measured 
scaling parameters A, B, and C are shown in Table IIVI 
As is to be expected in both cases the edge and cor- 
ner energies are more negative than the bulk bond en- 
ergies because the bonds at the edges and corners are 
less frustrated. The de-frustration effects can be seen to 
be ~ 50% stronger in the ±J case as compared to the 
system with Gaussian interactions. 

The Gaussian values can be compared with explicit 
bond-by-bond sample-by-sample measurements of bulk, 
edge, and corner bond energies. Averaging over 100 
L = 100 ground states with Gaussian bonds we obtain 
-0.6576(4), -0.698(2), and -0.66(2) for the mean local 
bulk, edge and corner bond energies, respectively. These 
values are compatible with the more precise values esti- 
mated by scaling, allowing for the fact that the definitions 
are identical for the bulk bonds but not for the edge or 
corner bonds. In these cases the scaling values implic- 



TABLE IV: Energy per bond for ff boundary conditions. 
The edge and corner energies are larger in magnitude than 
the bulk values because the bonds at the edges and corners 
are less frustrated. 



bond 


Gaussian 


±J 


Bulk 


-0.65740(1) 


-0.70099(1) 


Edge 


-0.7108(2) 


-0.7925(4) 


Corner 


-0.727(1) 


-0.821(4) 



itly include energy changes at slightly perturbed further 
bonds. 



V. PERIODIC-PERIODIC BOUNDARY 
CONDITIONS (PPBC) 

Periodic-periodic boundary conditions are the stan- 
dard geometry in which simulations are conventionally 
carried out, and to which the critical scaling rules dis- 
cussed in Sec. IIIII [Eq. ©] should apply. There are no 
edges or corners and all sites are on average in equiv- 
alent environments. The infinite-size limit energy per 
spin should be identical to the infinite size bulk energy 
per spin in the ffbc geometry. 

For the Gaussian case, high precision domain-wall 
measurements give excellent scaling with an exponent 
e = -0.282(2)fifliii and we should expect 



(19) 



up to correction terms. In fact we go further and test 
much more stringent assumptions: first, that corrections 
are negligible, and second that the energy per spin at 
large and moderate L extrapolates exactly to = at 
L = 1. Thus we test the relation 



l-L-(2-») 



(20) 



Here we have no free parameters as we adopt the value 
Coo = —1.31479 obtained above in the ffbc measurement, 
and the value 9 = —0.282 obtained from domain-wall 
measurements. We exclude i = 2 as this size is patho- 
logical in ppbc because of wrap-around effects for the 
interactions. For the "fit" = 0.92 for sizes from L = 3 
to 64. The data are represented in Fig.Owith a difference 
plot in the inset. It can be seen that this parameter-free 
expression represents the high precision data to within 
the statistical accuracy. Alternatively, assuming that the 
value is correct and that there are no corrections to 
scaling, we can fit leaving as a free parameter. We 
obtain a best fit for the stiffness exponent. 



e = -0.281(7) 



(21) 



i.e., v ~ 3.55, with x = 0.78. This is an ac- 
curate independent value in excellent agreement with 
the domain-wall estimates, as well as finite-temperature 
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FIG. 3: Energy for the ppbc Gaussian case. The sohd hne 
corresponds to a "fit" with no free parameters to Eq. I20II 
which yields a quahty of fit ~ 0.92. The inset shows the 
difference between the fit and the actual data. 
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FIG. 4: Energy for the ppbc ±J case. The solid line corre- 
sponds to a parameter-free "fit" according to Eq. 12UII with 
fioo = —1.401938 and ^ = 0. The inset shows the differences 
between the fit and the actual data. 



estimates li^iiSi The assumptions that we have made are 
either exact or are very good approximations. 

For the ±J case (Fig.0J the situation is less clear cut. 
If we make the same stringent parameter-free assump- 
tions, taking e^o = —1.401938 and = we obtain a 
rather mediocre fit, w 15 for the data from L = 3 
to L = 32 or w 6.5 for the data from L = 4 to 
L — 32. Relaxing the condition eL=i = gives a bet- 
ter fit, CL = eoo + 1.319/L2 with ~ 2.6 for L = 4 to 
L = 32. This could be understood as an analytic cor- 
rection term appearing. Such a term should be expected 
to also have the form L^^ and so in the ± J case where 
9 = Q it would simply modify the prefactor of the scaling 
term which is also proportional to . Alternatively, if 
we assume that for the range of L over which the energy 
differences can be measured the effective value of 9 is dif- 
ferent from zero, i.e., when incorporating the corrections 
to scaling into it, we can fit to Eq. (|20|l with as a free 
parameter. The best fit in this case from L = 4 to L = 32 
corresponds to = —0.04(1), with w 2.0. This discus- 
sion indicates that the strict zero free parameter scaling 
that gives virtually perfect agreement in the Gaussian 
case explains the form of to much lower precision in 
±J case. This could mean that unknown higher-order 
corrections to scaling must be taken into account to rep- 
resent the behavior also at the relatively small sizes ac- 
cessible for the ppbc case. We believe that this behavior 
is due to the discreteness of the bond distribution, lead- 
ing to a high ground-state degeneracy of the ± J model; 
see the discussion at the end of the next section. 



VI. PERIODIC-FREE BOUNDARY 
CONDITIONS (PFBC) 

We study a cylindrical geometry where the boundary 
conditions are free along one pair of edges, and periodic 
for the perpendicular pair. In the large-L limit the bulk 
site or bond energy term per spin or per bond should 
again be exactly the same as in the ffbc case. We as- 
sume that the free-edge size effect per edge site also has 
the same value as in the ffbc case. There are no cor- 
ners. Because the periodic boundary condition imposes 
a constraint, a priori one can also expect terms of the 
same type as in ppbc. The edge effects being allowed 
for, we assume that the scaling term is again propor- 
tional to L"*^^"^' and we allow for an analytical cor- 
rection term proportional to . The total number of 
bonds is 21? — L. Following the same arguments as in 
the ffbc case presented in Sec. II VI we test the expression 

eL = A* + ^{B*-A*)-^B*+C^L-^^-'^ + ^ , (22) 

where A* and B* are taken directly from the ffbc analysis, 
so we have only two free parameters for the fit. 

For the Gaussian pfbc data, the fit to points from L ~ 
3 to i = 320 is excellent, = 0.26 for Ci = 1.019(9) 
and £2 -0.203(7). The data are shown in Fig. EJ If 
there is no L^^ correction term, the fit is significantly 
worse. We note that for the pfbc geometry an analytic 
correction term is necessary, and that the fit does not 
extrapolate exactly to = at L = 1. The fit is thus 
less aesthetically pleasing than was the case for either 
ffbc or the ppbc cases; we can conclude that nevertheless 
a high quality consistent analysis can be made of the pfbc 
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FIG. 5: Energy for the pfbc Gaussian case and fit to Eq. H22|l 
(see the text). The inset shows the differences between the fit 
and the actual data. 
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FIG. 7: Domain-wall energies for Gaussian as well as ±J 
bonds. While in the Gaussian case the asymptotic power-law 
behavior appears for small system sizes, in the ±J case the 
asymptotic behavior where 6 = does not appear until lattice 
sizes L ~ 100. 
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FIG. 6: Energy for the pfbc ±J case and fit to Eq. I|22fl (see 
the text). The inset shows the differences between the fit and 
the actual data. 



Gaussian data over the entire range of L using physically 
reasonable assumptions for the fits. However, it is not 
clear to us why an analytic correction term is needed for 
the pfbc geometry while the ppbc fit is excellent without 
any such correction. 

For the ±J data, see Fig. El we proceed in exactly the 
same way. In this case the large- L stiffness exponenljiS 
6* = 0, so one should expect both the scaling and the 
analytical correction to contribute to a single composite 
term in L~^. With this assumption the fit to the pfbc 



data from L = 3 to L = 480 is very poor (x^ « 30). 

To understand this behavior, we can note that the 
domain-wall stiffness measurements— on the ±J ISG 
taken in the same pfbc conditions show strong and 
nonorthodox corrections to scaling which extend to L 
100; the data only attain the limiting behavior 6* = for 
very large L, as can be seen in Fig. [3 In the Gaussian 
case both the domain-wall stiffness and the size depen- 
dence of the energy in the pfbc geometry follow consis- 
tently the expected scaling rules. For the ± J case in pfbc 
geometry on the other hand corrections to scaling of the 
standard form are not adequate to explain the deviations 
from scaling which exist up to large values of L. This is 
probably due to the discreteness of the interactions for 
the ±J case. This discreteness leads to a high degen- 
eracy, which in turn allows in two dimensions to form 
domain-walls with almost no energy (as for the Gaus- 
sian case). But considerably larger system sizes have to 
be studied to find these low-energy paths, because, in 
contrast to the Gaussian case, there are no bonds which 
can be broken at the cost of very small energies. Every 
broken bond costs an energy 2 J. Indeed, by introducing 
possible zero bonds, i.e., treating a diluted system, the 
effect is reducediii As we will show in the next section, 
the scaling of the ground-state energy is related to the 
physical appearance of domain- walls, hence the same un- 
orthodox and strong corrections are to be expected for 
the ± J case. 
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FIG. 8: DifTerence between the ground states of one realiza- 
tion (L — 100) for free and periodic boundary conditions in 
the a;-direction. Bonds are changed along the vertical axes. 



VII. APPEARANCE OF DOMAIN WALLS 



TABLE V: Number of samples for the ground-state calcu- 
lations for the samples (Gaussian distribution of the bonds) 
having first full free, then free boundary conditions in the 
y-direction and mixed random- weak/inverted-hard boundary 
conditions in the x-direction (with a fraction p of the bound- 
ary bonds being inverted hard), for different system sizes L 
and values of p. For the difi'erent values of p > 0, always the 
same number of samples are taken for a fixed size, except for 
the largest size L = 160. 



_L p^O p > 

6 21000 1000 

10 21000 1000 

16 20000 

20 10000 1000 

30 20000 

40 20000 1000 

60 15000 

80 13000 1000 

120 15000 

160 15250 700-4600 



We have seen, in particular for the model with a Gaus- 
sian distribution of the interactions, that for free bound- 
ary conditions the size dependence of the ground-state 
energy can be explained by simple algebraic terms. As 
soon as periodic boundary conditions are applied, a term 
^ appears, 9 being the exponent describing the 

size-dependence of the energy of elementary excitations 
such as droplets or domain walls. In this section we want 
to therefore test whether really the occurrence of domain 
walls, induced by the boundary conditions, is respon- 
sible for this scaling. For this purpose, we study sys- 
tems with Gaussian interactions and ffbc, calculate the 
ground state, then we change the boundary conditions to 
pfbc (by adding one column of random bonds wrapping 
around the system), recalculate the ground state, and 
compare the two different ground states. In Fig.|Slthe re- 
sult for one sample of size L = 100 is shown. The result is 
typical: no system-spanning domain walls are obtained, 
but a collection of smaller and few larger droplet-like ex- 
citations pinned to the boundary. We observe that quite 
often the larger droplets span even larger fractions of the 
systems. 

To test whether system-spanning domain walls play a 
role in this framework, we perform the following simula- 
tion: We start with ffbc, calculate the ground state, and 
then switch to pfbc. Now a fraction p of the bonds which 
wrap around the system is chosen as inverted hard bonds. 
This means that they have a very large magnitude, such 
that they are satisfied in any ground state. The sign is 
chosen in a way such that it forces the two adjacent spins 
to take opposite relative orientations with respect to the 
ffbc ground state, i.e. exactly one spin flips. This means, 
the inverted hard bonds force a domain wall into the sys- 
tem at this position. The remaining fraction (l — p) of 
the additional bonds is again chosen from a Gaussian dis- 
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FIG. 9: (Color online) Probability Pdw of the occurrence of 
a system-spanning domain wall as a function of the fraction 
p of inverted hard bonds wrapping around the system, for 
different system sizes L. The lines connect the dots and are 
guides to the eye only. The inset shows Pdw(p = 0) as a 
function of L. We see that for p = in average only ~ 10% 
of the domain walls are system spanning. 



tribution with zero mean and standard deviation unity. 
The number of samples studied as a function of L and p 
are shown in Table IVl 

In Fig.|51the probability Pdw that a domain wall spans 
from the top to the bottom is shown as a function of p. 
Clearly, Pdw increases with growing p and reaches unity 
for largest system sizes when p — > 1, as expected. For 
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FIG. 10: (Color online) Length fow (in units of L along the y- 
direction) of the longest nontrivial domain wall as a function 
of the fraction p of inverted hard bonds wrapping around the 
system, for different system sizes L. Lines are guides to the 
eyes only. The inset shows luw{p = 0) as a function of L. 

large p, Pdw grows with system size. Near p = 0, no such 
trend is visible, see also the the inset of Fig. |5] In par- 
ticular, no clear crossing of the curves can be observed, 
which would be an indication that a finite fraction of in- 
verted hard bonds is necessary to create system spanning 
domain walls. Thus, system-spanning domain walls exist 
in the thermodynamic limit also for p = 0, but not many. 

To obtain a clearer picture, we also measure for each 
system the size of the largest nontrivial domain wall (the 
largest trivial domain contains the spins which have not 
changed) in units of the system size L, when projected 
onto the y axis. The resulting length Zdw averaged over 
disorder is displayed in Fig. ^| On can see that the 
domain walls grow with increasing fraction p of the in- 
verted hard bonds. For p = (see the inset of Fig. [TU^ . 
?DW seems to converge to a quite large value, larger than 
0.5. This means that the periodic boundary conditions 
create domain walls of order system size, although not 
spanning the system. For these domain walls (one can 
call them also droplets pinned at the border) indeed a 
scaling of their energy proportional to can be ex- 
pected. Hence, there is a intuitively clear reason for the 
appearance of the nontrivial scaling term, which we have 
found above when studying the ground-state energy of 
two-dimensional systems. 

VIII. CONCLUSIONS 

We have analyzed the size-dependent energies for two- 
dimensional ISGs with Gaussian and ± J interactions, up 
to large sizes and for three different types of boundary 



conditions: free along both axes of the sample (ffbc), 
periodic (ppbc) and free/periodic (pfbc). 

With a simple scaling expression and three free pa- 
rameters, excellent fits can be obtained to the ffbc free 
boundary condition data, from the isolated spin L — Inp 
to large sizes, L = 320 for the Gaussian and L — 480 for 
the ±J cases. The fitting parameters can be translated 
into bulk-, edge-, and corner-spin (or bond) energies. The 
bulk ground state energy per spin is Coo = —1.31479(2) 
for the Gaussian and Coo = —1.40197(1) for the ±J case. 
The former is considerably more accurate than previ- 
ous estimates, and the latter is consistent with the high- 
precision published value The effective edge and corner 
bond energies do not seem to have been measured before; 
they are more negative than the bulk values as might be 
expected because of geometrically reduced frustration. It 
is remarkable that although the entire range of L includ- 
ing L = 1 is used, three parameters only are sufficient for 
an excellent fit in both cases. 

For the periodic-periodic boundary condition geome- 
try, a scaling analysis of the Gaussian system data can 
be made using only parameters (coo and the stiffness ex- 
ponent 0) whose values are already accurately known. 
The zero free-parameter "fit" passes through eL=i = 
just as for the ffbc case. Thus for the Gaussian system 
the scaling in these two geometries appears to be excel- 
lent, implying that any correction terms present are very 
small. 

For the ± J system, the ppbc fit is less satisfying than 
for the Gaussian case, as an acceptable fit requires the 
introduction of a correction to scaling term. Even with 
this supplementary term included the global fit is not as 
good as for the Gaussian case. 

In the mixed pfbc geometry the fit to the Gaussian 
data needs the introduction of a weak analytic correction 
term, and the extrapolation to L = 1 does not quite 
go through eL=i = 0. The scaling is thus sHghtly less 
aesthetically pleasing than for the ffbc or ppbc geometries 
but follows standard correction to scaling rules. 

For the ± J system any fit to the pfbc data based on the 
same approach as for the Gaussian is very poor. In this 
geometry, for the size dependence of the energy as for the 
directly measured domain-wall stiffness, there are strong 
deviations from scaling that do not follow the orthodox 
correction to scaling behavior. As these deviations do not 
appear in the Gaussian case, we associate them with the 
degenerate ground state of the 2c? ± J ISG, as discussed 
above. This hypothesis could be checked by measure- 
ments on other 2d ISG systems with degenerate ground 
states. 

By studying explicitly sample by sample the same sys- 
tem having ffbc and pfbc, we have seen that indeed do- 
main walls of order of the system size are created by the 
periodic boundary conditions, although usually these do- 
main walls (or pinned droplets) do not span the whole 
system. The appearance of these domain walls explains 
the occurrence of the term in the size dependence of 
the ground-state energy for periodic boundary conditions 
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in a quite natural way. 
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